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We demonstrate the first implementation of a quantum algorithm on a liquid state nuclear mag- 
netic resonance (NMR) quantum computer using almost pure states. This was achieved using a 
two qubit device where the initial state is an almost pure singlet nuclear spin state of a pair of ^H 
nuclei arising from a chemical reaction involving para-hydrogen. We have implemented Deutsch's 
algorithm for distinguishing between constant and balanced functions with a single query. 
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INTRODUCTION 



The practical implementation of quantum algorithms 
has so far been dominated byliguid state nuclear mag- 
netic resonance (NMR) H E Q techniques H B S 
0, S IE U^- The many algorithms successfully im- 
plemented on liquid state NMR systems include, most 
notably, Deutsch's algorithm [Q, Qi Grover's quantum 
search p,lg|, and Shor's factoring algorithm jl3|. This 
unrivalled success can be attributed to the high degree of 
coherent control over small spin systems and also to the 
long decoherence times in comparison with the timescales 
of the computation. However, one of the problems haunt- 
ing the prospects of liquid state NMR as a practically 
scalable and "truly quantum" implementation, was that 
the states encountered in almost all of the previous im- 
plementations were highly mixed jT^, [l3 • 

The highly mixed states used in NMR raised two ma- 
jor concerns. First, there was the issue of initializing the 
spins into a well-defined state [ij, [lj| • It was shown that 
one could circumvent this problem [3) 13 by preparing 
pseudopure states [l^ [la, [l3 , which behave effectively 
(up to a certain scaling factor) as pure states. This ap- 
proach worked well for small spin systems, but cannot 
be extended to larger systems, because of the exponential 
loss in signal with increase in the number of spins Il2lll4| . 
Secondly, the highly mixed states encountered in conven- 
tional liquid state NMR are separable ,11]. As entan- 
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glement could not exist in these systems, some authors 
suggested that NMR quantum computing could perhaps 
be explained using classical models llSj. With separa- 
ble states, even the quantum nature of NMR quantum 
computing was under question! 

Motivated by these problems, we have previously 
demonstrated [I3 the preparation of an almost pure two- 
spin state, which lies above the entanglement threshold 
|23,'2J| , for the two hydride "'^H nuclei in the organometal- 
lic compound Ru(H)2(CO)2(dppe), where dppe indicates 
l,2-bis(diphenylphosphino)ethane. This was achieved 
using laser induced addition of pure para-hydrogen to 
Ru(C0)2(dppe) as discussed below. We now describe an 
implementation of the Deutsch algorithm on this spin 
system, demonstrating that coherent manipulations of 
spins in our pure spin system can also be carried out. 

Our experiments solve the two-fold problem of initial- 
ization and separability. We bypass the need for prepar- 
ing pseudopure states, as our computation starts off di- 
rectly with an almost pure Werner state .22j . Further- 
more, the initial purity of the system lies above the en- 
tanglement threshold. We do not claim that our imple- 
mentation actually involves entangled states — only that 
entanglement can be distilled from these states. The use 
of entanglement proper in the Deutsch's algorithm has 
been discussed in more detail elsewhere 23]. In short, 
the current work constitutes the first demonstration of a 
quantum algorithm in liquid state NMR using an almost 
pure initial state. 

The remainder of the paper is organized as follows. In 
section ^ we describe the use of para-hydrogen to gen- 
erate spin systems in almost pure initial states, while in 
section llTT] we give a brief summary of Deutsch's problem 
and the quantum algorithm devised to solve it. Liquid 
state NMR can be used to implement the algorithm as 
with our two spin pure state; this is discussed in section 



IrVl and the experimental procedures and results are de- 
scribed in section [3 Finally, we present a concluding 
summary in section IVII and also outline some possible 
directions for future work. 
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II. P^flA-HYDROGEN 



FIG. 1: Disentangling network, rotating Bell states into eigen- 
states. 



The high spin-state purity in our experiment is a re- 
sult of an effect called para-hydrogen induced polariza- 
tion (PHIP) in I2I El Ell- The existence of the para 
spin isomer of dihydrogen molecules, H2 , is a consequence 
of the Pauli principle and the symmetrization postulate 
128]. These require the overall wavefunction of the H2 
molecule to be antisymmetric with respect to interchange 
of the fermionic ^H nuclei. The translational and vibra- 
tional wavefunctions of H2 are always symmetric, as is 
the ground state electronic wavefunction, and so over- 
all antisymmetry is achieved by choosing the rotational 
and nuclear spin wavefunctions such that their product 
is antisymmetric. It follows that H2 molecules in even 
rotational states (J = 0, 2, . . .) possess an antisymmet- 
ric nuclear spin wave function and correspond to nuclear 
spin singlets, termed para; similarly, H2 molecules in odd 
rotational states (J = 1, 3, . . .) possess a symmetric nu- 
clear spin wave function and correspond to nuclear spin 
triplets, termed ortho. 

Interconversion of the two spin-state isomers of hydro- 
gen is normally forbidden by angular momentum selec- 
tion rules, but in the presence of a paramagnetic catalyst 
the symmetry of the H2 system is broken, allowing these 
rules to be overcome. Thus if hydrogen is cooled in the 
presence of an appropriate catalyst, it will become en- 
riched in the (low energy) para form. Upon moving away 
from the catalyst interconversion is again suppressed, and 
so the ortho / para ratio in effect remembers the temper- 
ature of the last conversion surface encountered. A tem- 
perature of 20 K is sufficiently low to cool 99.8% of H2 
molecules into the J = state |26j and hence produce 
essentially pure para-hydrogen. 

The sole antisymmetric state of two -'^H nuclei is the 
singlet state, and so pure para-hydrogen will also have 
a pure nuclear spin singlet state. This ability to pre- 
pare a pure initial state is obviously attractive for NMR 
quantum computing experiments, but the para-hydrogen 
molecule cannot be used directly for NMR quantum com- 
puting, as it is NMR silent due to its high symmetry. 
This can be overcome by means of a chemical reaction, 
producing a new molecule, in which the two hydrogen 
atoms can be made distinct (I and S) and can be sepa- 
rately addressed. In conventional para-hydrogen experi- 
ments the chemical reaction is slow in comparison with 
the frequency difference between the I and S spins of 
the reaction product. This causes the off-diagonal terms 
in the density matrix to dephase, producing a separable 
state. This dephasing can be overcome (neglecting re- 
laxation) |29j | by applying an isotropic mixing sequence 
[30|, but a much simpler approach is to ensure that the 



reaction is rapid in comparison with the dephasing and 
relaxation timescales. We achieve this through addition 
of para- hydrogen to a highly reactive species, generated 
by laser flash photolysis of a stable precursor. 

The state that we actually prepare in our para- 
hydrogen experiments is observed to be a Werner singlet 
state ^21 of the form 
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where \4> ) is the singlet vector, which can be written in 
the computational basis as 



|V'-)= ^(101) -|io)). 



(2) 



and the polarization e determines the purity, with e = I 
corresponding to the pure state; for our system [l9j e « 
0.92, and so the state is almost pure. The singlet state is 
one of the four Bell states [3l|, |(/)±) = (|00) ± |ll))/%/2 
and IV'^) = (|01)±|10))/-\/2, all of which can be intercon- 
verted using local unitary transformations |8ll l83 | . The 
Bell states can also be converted to basis vectors in the 
computational basis (the usual starting point for quan- 
tum computations) using a network of the form given in 
Figure ^ which is a disentangling network. The output 
qubit I a) is determined by the label {(f) or ip) of the input 
state lipin) and \b) is determined by the sign (-1- or — ); 
eigenstates can be converted back to the Bell states, us- 
ing the reverse of the circuit in Fig. ^ and the mapping 
can be compactly expressed as 
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where pb) is the classical NOT of \b). With circuits 
built around the disentangling circuit, we can achieve 
conventional two qubit initial states, starting off from 

\ll^Z7l) = IV'")- 



III. THE DEUTSCH PROBLEM 

Deutsch's problem [33 was the first problem to be 
solved by a quantum algorithm. Although its useful- 
ness as a real-life problem is limited, it is nevertheless 
a convenient testing ground for quantum computation. 
A refined version of Deutsch's problem [33| is described 
in the following way. 
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TABLE I: Truth table for the four single-bit functions. 
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Suppose we are given a binary string s — {0, 1}" of 
length n and a binary function / acting on s mapping it 
onto a single bit, or 1, i.e., f : {0, 1}" -^ {0, 1}. Fur- 
thermore, we are promised that / is either constant or 
balanced. (A constant / outputs the same value, either 
1 or 0, for all 2" possible input strings s, whereas a bal- 
anced / outputs for exactly half of the strings and 1 
for the other half.) The problem is then to find whether 
/ is constant or balanced using the minimum number of 
queries. In the classical setting, answering this question, 
in the worst case, would require (2"/2) + 1 queries. How- 
ever, with a quantum algorithm it is possible to solve this 
problem in a single query, giving an exponential improve- 
ment over the worst-case classical version. 

The simplest form of Deutsch's problem concerns the 
case of n = 1, with functions mapping a single bit onto 
a single bit, / : {0, 1} -^ {0, 1}. There are four such 
functions, their operation being described by the truth 
table given in Tabled Out of these four functions, two 
are constant (/oo and /n) and two are balanced (/oi and 
/lo). Classically the problem is solved by evaluating first 
/(O) and then /(I) and comparing the results. Equiva- 
lently, / is constant if /(O) © /(I) = and balanced if 
/(O) © /(I) = 1, where © indicates addition modulo 2. 
Deutsch's algorithm enables us to find this global prop- 
erty of / in a single query. 

The conventional approach to solve the single-bit prob- 
lem is to use two input and two output qubits and then 
implement / using carefully selected reversible gates. 
These gates are represented by the propagators C//, 
whose action on the input qubits can be expressed as 



FIG. 2: Quantum circuit for the classical evaluation of f{x) 
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leaving the first qubit unchanged. If |y) is initialized as 
|0), the second qubit holds the value of f{x) after the 
computation, 



j:)\0)^\x)\f{x)). 



(5) 



There are four propagators, Uqo, Uqi, Uiq and Uu, cor- 
responding to the four different functions. A suitable 
choice of propagators is 



C/oo = 1 «) 1 


(6) 


c/oi = |o)(o|(»a + |i>(i|®(7. 


(7) 


t/io = |0)(0|(^fT, + |1)(1|®]1 


(8) 


C/ii = 1 (g) <j^ 


(9) 



where 1 is the unit matrix of order 2 and cr^ is a Pauli 
spin matrix [3l| . It is straightforward to check that with 
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FIG. 3: Quantum circuit for solving Deutsch's problem. H 
represents a Hadamard gate. 



appropriate initial states these propagators indeed repro- 
duce the transformation given in Equation|31 In the usual 
language of quantum networks C/qo is a do-nothing gate; 
Uoi is a l-controUed-NOT gate; Uio is a 0-controlled- 
NOT gate and Uu is a NOT gate on the second qubit. 
(Reference [35j contains a concise description of quantum 
logic gates in NMR). 

We can build a quantum circuit from these propaga- 
tors, which is depicted in block form in Figure |2| With 
this network, the second qubit evaluates f{x), although 
the evaluation is classical and we still need two separate 
evaluations to determine the nature of /, once with x = 
and once with x — 1. This circuit is, however, useful 
for illustrative purposes. For example, using this circuit 
one could verify the truth table in Table ^ To solve 
Deutsch's problem in a single step, however, we exploit 
superposition states, preparing them from the eigenstates 
using Hadamard gates [31|. The corresponding network 
is given in Figure O and the algorithm proceeds as fol- 
lows: the computation starts off in the state |0) (g) |1), 
the Hadamard gates prepare the superposition state, fol- 
lowed by the unitary transformation Uf, and the final 
Hadamard gates then convert the superpositions back to 
eigenstates, with the first qubit ending in /(0)©/(l) — 
the solution to the problem! 



IV. IMPLEMENTING THE ALGORITHM 

We have implemented the Deutsch algorithm on an 
almost pure (e « 0.92) two qubit state using an NMR 
quantum computer. We also evaluated f{x) classically 
using the procedure already outlined. In an NMR quan- 
tum computer, the qubits are spin-1/2 atomic nuclei (in 
this case ^H nuclei) placed in a strong magnetic field, 
and the gates are implemented through a sequence of 
radio- frequency (RF) pulses, interspersed with properly 



timed delays, during which the spin system evolves un- 
der its own Hamiltonian. Details may be found in reviews 
of quan tum computing with liquid state NMR, such as 
I35, 36, j7, ^8]. Here we use product operator notation 
[i,^, 40] for the states and propagators. Pulse sequences 
are written from left to right and we label the two spins 
as / and S, in keeping with NMR parlance. The notation 
9a is used to represent a rotation by 9° about the a axis. 

To implement the classical and quantum versions of 
the problem (Figures |21 and |2Jl, the individual gates have 
to be built from RF pulses and frame rotations. For ex- 
ample, the Hadamard gate on a single qubit is a 180° 
rotation about a tilted axis |35j and can be achieved us- 
ing the sequence 180z 90-y. Pulse sequences can also be 
simplified based on the properties of product operators: 
for example, if the Hadamard acts on an eigenstate, the 
initial z-pulse can be dropped out as these states com- 
mute with z-pulses. We are therefore simply left with 
90-y pseudo-Hadamard operations which can be used as 
a simpler replacement for proper Hadamard gates. Like- 
wise, if a computation ends in a diagonal state, we can 
shift the 180^ pulses to the end (right) of the sequence, 
changing, of course, the phase of the 90° pulses, and col- 
lapse the z-pulse with the final state. We are then left 
with a simple 90y pulse in place of a strict Hadamard. 
In the same spirit, the z-rotations can be moved within 
a sequence, either to the left or to the right, and can be 
viewed as rotations of the reference frame, often referred 
to as the abstract reference frame IT7|. We have used 
the abstract frames approach together with composite 
2-rotations HJ] to simplify our sequences. 

The unitary transformations corresponding to the four 
possible t//'s are given in Equations 1^1 f/00 is an iden- 
tity gate, and translates to "do nothing", while Uu is 
a selective NOT gate on the second qubit, 1805j;. The 
remaining two transformations, t/01 and Uio can be ex- 
pressed in product operator notation as, for example. 



a b 


Pulse sequence 





A = [Ti]90y [t2] 180. [ra] 180,, [ri]90. 


1 


B= [Ti]90y[T2]1804r2]90y 


1 


C= [ri]90^[r2]180.[r2]90-y 



and 



90S-y 90/^ 905*^ (-90)2/^5'^ 90Sy (10) 



90S-y 90I-Z 90Sz {90)2hSz 90Sy (11) 



respectively. These transformations involve frame rota- 
tions, evolution under scalar coupling terms (2IzSz) and 
soft pulses 3J, selectively exciting the S spin. Selective 
pulses were implemented using the Jump and Return 
pulses described previously |42l |. so that our final sim- 
plified sequences comprise only hard pulses and delays. 
The sequences for the non-trivial transformations, Uqi, 
Uio and Un are 

Poi = 90_,[ti]9045[2ti + T2]180,[T2]90i35[ri]90_, (12) 
Pio = 90,[ri]9045[2Ti + T2]180,[r2]90i35[ri]90, (13) 
Pii=90j,[2ri]90_j,90, (14) 

where ri = 1/46 and T2 — 1/4J, S being the frequency 
separation between the spin resonances and J the spin- 
spin coupling constant, both measured in Hz. 



TABLE II: Pulse sequences for mapping the singlet state \^p ) 
onto the input states \ab). 



We have so far overlooked another important issue, the 
preparation of initial states. The para-hydrogen derived 
state is an almost pure singlet, but the algorithms expect 
the quantum computer to begin in some eigenstate \ab), 
with the values of a and b depending on the algorithm. 
With circuits built around the disentangling circuit. Fig. 
n we can achieve the desired input states, starting off 
from |-0i„) = \'ip~). The simplified pulse sequences to 
accomplish this transformation are shown in Table UTI 

Implementation of the classical and Deutsch algo- 
rithms involves preparing the almost pure state |'0~)j 
converting it into the desired initial state using a recipe 
from Table m and applying the appropriate pulse se- 
quence for the algorithm. Finally, we need to read out 
the result at the end of the computation. At this stage 
the output states are diagonal, and so cannot be directly 
observed, but a simple hard 90^ acquire pulse rotates 
these states to the measurement basis: the |0) state gives 
an NMR signal I^ (positive absorption) and |1) gives a 
signal —Lj: (negative absorption). The output states can, 
therefore, be unambiguously assigned as |0) or |1) by ex- 
amining the two multiplets. 



V. THE EXPERIMENT 

The two qubit system for our NMR quantum computer 
comprises the two hydride ^H nuclei in the organometal- 
lic compound Ru(H)2(CO)2(dppe), where dppe indicates 
l,2-bis(diphenylphosphino)ethane and the hydride hy- 
drogen atoms are derived from para-hydrogen. We pre- 
pared essentially pure para-hydrogen at a temperature of 
18K in the presence of a charcoal-based catalyst. The 
gas was introduced into a 5 mm NMR tube contain- 
ing the precursor compound Ru(C0)3(dppe), dissolved 
in dg-benzene. The preparation of the precursor and 
finer experimental details are identical to previous work 
[Slillil- The NMR tube was then transferred into a 
400 MHz spectrometer fitted with a -'^H/'^^P probe mod- 
ified for in situ photolysis |44| . The spectrometer trig- 
gered an MPB Technologies MSX-250 pulsed XeCl ex- 
cimer laser which fired a 12 ns UV pulse of wavelength 
308 nm, irradiating the active region of the NMR sam- 
ple and producing the unstable species Ru(C0)2(dppe). 
This unstable intermediate reacts with dissolved para- 
hydrogen on the sub-microsecond timescale |45l | leading 
to the product of interest, Ru(H)2(CO)2(dppe). The two 
hydrogen nuclei which inherit the pure singlet spin state 



Experiment 



Pulse sequence 



Classical /(O) A—G—Uf—90y 

Classical /(I) C— G— ?7/— 90^ 

Quantum /(O) e /(I) B— G— 90^— ^7/ 

TABLE III: Pulse sequences for the Deutsch algorithm. Gra- 
dient fields are represented by G. 



from the para-hydrogen comprise the two qubits of our 
quantum computer. 

The hydride resonances appear at — 7.55ppm (spin /) 
and — 6.32ppni (spin S), with a frequency separation of 
d — 492 Hz. The ^H transmitter frequency was placed 
exactly between the two resonance frequencies, and cou- 
plings to "^^P nuclei were removed by GARP decoupling 
[4|5j applied continuously throughout the experiments. 
The Ti and T2 constants for the hydride peaks were mea- 
sured to be 1.7 and 0.58 s respectively while the hydride J 
coupling (■^ Jhh) was 4.6 Hz. The laser flash acts as an ini- 
tialisation switch, generating the pure state IV'in) = |V'~) 
on demand, which is subsequently used for the implemen- 
tation. The pulse sequences for the classical and quan- 
tum evaluations are described below and summarized in 
Table iHll 

Using our system, we were able to analyze /(O) and 
/(I), corresponding to the classical evaluation of these 
functions (see Figure [2)). For example, for the determi- 
nation of /(O) the initial state |00) can be prepared from 
the singlet state \ip~) using the sequence A (see Table 
Hlj) . At this stage, we clean our state by applying a strong 
field gradient |3, llfl ; this will not affect the component 
of the state in the desired eigenstate, but the majority of 
unwanted terms will be dephased Q- We then apply the 
transformation Uf and finally measure with a 90j, acquire 
pulse. The resulting spectra for the four possible choices 
of Uf are shown in Figure 0| From Table U we observe 
that /(O) = for /oo and /oi and /(O) = 1 for /lo and 
/ii; the result is encoded in the state of spin S (on the 
left), whereas spin / (on the right) remains in the state 
|0). Thus we should obtain two kinds of spectra, cor- 
responding to Ix i Sx, with the / resonance in positive 
absorption and the S resonance in positive or negative 
absorption, depending on the value of f{x). The two 
kinds of spectra are clearly seen (Figure 2J). 

In the same way we evaluated /(I) by converting the 
singlet state into |10) using the sequence C and a strong 
gradient field, applying the operation Uj and the acquire 
90y pulse, and observing the result; spectra are shown in 
Figure [5| We again observe two kinds of spectra, — I^: ± 
Sx '■ spin / is always in negative absorption and S changes 
phase, encoding the value of /(I). Our results are in 
complete accord with Table Q] 

Finally, we implemented the quantum Deutsch algo- 
rithm using the sequence B, followed by the gradient, a 
pseudo-Hadamard and the transform Uf. An additional 
simplification occurs in the quantum version: the final 
pseudo-Hadamard can be cancelled with a 90y acquire 
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FIG. 4: Results of the algorithm for the classical determina- 
tion of /(O). (a) The result of applying Uoo- The spectrum 
corresponds to I^c + Sx and both spins are in positive ab- 
sorption. This spectrum can also serve as reference for the 
phasing, (b) The result of applying f/01 , corresponding again 
to Ix + Sx- (c): The result of applying Uiq. The spectrum 
now corresponds to Ix — Sx, the left spin / is still in posi- 
tive absorption and the right spin S is now inverted, (d) The 
result of applying Uu, corresponding to Ix — Sx- Spin / al- 
ways remains in positive absorption whereas S changes sign 
according to /(O). 
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FIG. 5: Results of the algorithm for the classical determina- 
tion of /(I). Labelling and interpretation are the same as 
in Figure 0] Spin / always remains in negative absorption 
whereas /(I) determines the sign of the signal from spin 5*. 



pulse and the two pulses can therefore be dropped al- 
together. The system is, therefore, already in the mea- 
surement basis after the transformation Uf and can be 
observed directly. The resulting spectra are shown in 
Figure El The observed states after the quantum evalua- 
tion are Ix — Sx and —Ix — Sx, for the constant (/ooi /ii) 
and balanced (/oi, /lo) functions respectively. 

Although the overall results are clear, our spectra show 
several imperfections. Each multiplet should in principle 
show a perfect absorption lineshape (either positive or 
negative), and all the multiplets in all the spectra should 
ideally have the same intensity. In fact dispersive com- 
ponents are clearly visible, and there are considerable 
variations in multiplet intensities, both within spectra 
and between spectra. These imperfections are most pro- 
nounced in the spectra involving f/01 and C/10, which were 
obtained with the longest and most complicated pulse 
sequences. The distortions can arise from many factors 
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FIG. 6: Results of the algorithm for the quantum determi- 
nation of /(O) © /(!)• Labelling and interpretation are the 
same as in Figure 0] Spin S always remains in negative ab- 
sorption while the sign of the signal from spin / encodes the 
final result. 



such as finite coupling evolution during delays, RF inho- 
mogeneities, miscalibrated pulses and decoherence. 

To simulate the effect of decoherence under both Ti 
and T2 relaxation processes, we employed a model based 
on the operator sum representation 10, "31] for phase and 
generalized amplitude damping. The model is crude and 
assumes independent uncorrelated relaxation processes, 
but is nonetheless useful in obtaining a rough first es- 
timate of the effects of decoherence. We observe that 
the damping processes should result not only in an over- 
all decrease in the signal intensity, but for the case of 
the Uqi and Uio propagators can also lead to intensity 
imbalances within a spectrum: for experiments involv- 
ing these two propagators one of the spins spends more 
time in the transverse plane than the other, and is there- 
fore more prone to phase damping. In the worst cases, 
this can cause an overall intensity imbalance of about 
2 : 3, similar to the observed ratios. However for a small 
two qubit system such as ours, the effects of decoherence 
are not big enough to obscure the measurement results. 
Comparing the intensities of the final spectra with the in- 
tensity of spectra obtained from the initial singlet state. 



we estimate the purity of the final state to be about 50 
to 60 percent, indicating that we remain above the en- 
tanglement threshold throughout the implementation. 



VI. CONCLUSIONS AND FUTURE WORK 

We have demonstrated an implementation of the 
Deutsch algorithm to distinguish between constant and 
balanced functions using an NMR quantum computer 
with an almost pure initial state. We have also used 
the same system for classical evaluations of the function 
values. Unlike previous implementations built around 
pseudopure states, our implementation uses almost pure 
states which can be made on demand by laser flash pho- 
tolysis. Our earlier work 19] suggested that we could get 
around the problem of low initial polarizations in liquid 
state NMR by making high purity states using the para- 
hydrogen approach. The present experiment is one step 
forward in showing that the high polarizations can be put 
to use in the form of a demonstrable quantum algorithm. 

It seems that NMR quantum computing with pure 
states could in principle be scalable, although more work 
is needed to verify the validity (or otherwise) of this 
claim. We are seeking to extend this work to larger spin 
systems and more complex algorithms using a combi- 
nation of direct para-hydrogen induced polarization and 
the transfer of spin-state purity between different spins. 
We are also investigating the possibility of using para- 
hydrogen techniques to reset qubits midway through a 
quantum computation, thus permitting repeated error 
correction schemes. 
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